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Abstract: In this paper we consider matrix and vector models in the large A^ limit (N x A' 
matrices and vectors with A^ 2 components). For the case of zero-dimensional model (D=0) 
it is proved that in the strong coupling limit g — > oo statistical sums of both models coincide 
up to a coefficient. This is also true for D = 1. 



1. Introduction 

Matrix models have applications in various branches of physics, but for us the main 
motivation will be the possibility of using them in quantum field theory. In particular, 
as it has been shown by 't Hooft [1] , with a proper normalization of coupling constants 
the expansion of Green's functions (and the partition function) in corresponds in 
terms of ribbon Feynman diagrams to the expansion in the genuses of Riemannian 
surfaces, on which this graph can be drawn without self-intersections. However, 
even if the genus of the surface is fixed, there are infinitely many corresponding 
ribbon graphs and in the majority of interesting cases, unfortunately, there's still no 
algorithm for their summation. 

The limit iV — > oo is also interesting from the point of view of the Ads/CFT 
correspondence [2]. In this case one should consider the large g limit of the sum of 
planar diagrams (N — > oo). This problem is far from being solved not only for gauge 
theories but also for the simpler case of scalar matrix models. Exceptions are cases 
D = and D — 1, where the algorithm of planar diagram summation has been found 
[3] for real and hermitian matrix models. Later (at D — 0) it was generalized for 
complex matrix models [4]. A review of the state-of-the-art condition of the theory 
of zero- dimensional matrix models can be found in [5]. 

Vector and matrix models at D — 0, 1 reveal remarkable similarity in the limit 
g — > oo. The dependence of the partition function on the coupling constant in this 
limit is the same for vector and matrix models. This suggests that such coincidence 
can also take place in higher dimensions. 

However, in this paper we will restrict ourselves to the consideration of a D = 
complex matrix model, i.e. we consider the following partition function: 

/I g 
dp dtp exp[--tr(v?V) - — tr((y>V) 2 )], 

where ip is a iV x N matrix. The notation dip' dip should be interpreted in the 
following way: 

dip" 1 dip = Dip = Y\ dRetpij dlmipij. 
Free energy is defined by 

_ lnZ BV 

We will show that at D = in the strong coupling limit (g — > oo) the leading asymp- 
totics of the free energy of this matrix model coincides with the leading asymptotics 
of the free energy of the corresponding vector model. At D — 1 this is also true, 
but a more thorough consideration of this case is necessary. Of course, especially 
interesting are the cases D > 2, but there the question of whether the described 
hypothesis is true is still open. We would like to emphasize from the beginning that 
at D > 2 ultraviolet divergences appear both in matrix and vector models, therefore 
it only makes sense to speak about properly regularized theories. 

Notice that in all dimensions linear vector models with interaction 
N 2 F \-±2 tr(</^ <p)) are exactly soluble, i.e. Green's functions of singlet variables can 
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be found in explicit form, the analogy between matrix and vector models which 
supposedly takes place can help in the analysis of (much more complicated) matrix 
models. 



2. D=0 bivector models 

For the comparison of matrix and vector models we first calculate the partition 
function of the vector model. Let's explain what we mean by vector and bivector 
models. The distribution p of random variables x a (a = 1...M) in a vector model, 
by definition, depends only on the square of the vector x: p = p{xl). 
Consider the partition function 



Z 



BV 



Dip exp 



■itr(^V) + iV 2 j p ( _l tr( ^ )) 



(2.1) 



Although formally this model is a matrix model, in fact it is a vector model, because 
if we introduce a vector y a with M = 2N 2 components 

(yi = Reip {il ,y N 2 +i = lm<p itl ,i < N), (y, t = Reipi„ Ni2 ,y N 2 +i = lm<pi- N j,N < i < 2N)..., 

then the distribution turns out to be just a function of y\. We call such models 
bivector models. They are characterized by the fact that their symmetry group is 
much wider than in the corresponding matrix model: for example, in this case the 
matrix model has a symmetry group U(l) x SU(N)l x SU(N)r (multiplication by 
a phase factor, left and right matrix multiplication x ), whereas the vector model is 
invariant under U(2N 2 ). This large symmetry is what lets us find an exact solution 
for (bi)vector models. 

First we want to find the N — > oo limit of the free energy. Inserting a unity into 
the integrand of (|2.1|) we obtain: 



Z 



BV 



Dip dp 5 (p — tr (<p'(p)) exp 



[2.2) 



Let us write the Fourier transformation for the delta-function: 



'BV 



,2,)-./ 



Dip dp dX exp 



-itr(pV) + N 2 F(j^p) + i\{p - tr(pV)) 



2.3) 

The integral over ip is gaussian, moreover it degenerates into a product of independent 
integrals because the components of the matrix ip are not linked: iviip^ip) = ^ |^jj| 2 - 



As a result, the integral over ip can be taken, and we get: 

iV 2 



'BV 



(27T) 



dpdX 



71 



1/2 + a 



exp 



N 2 F{—p)+i\p 



(2.4) 



1 Multiplication by a phase factor does not reduce to a combination of left and right multiplication by 
unimodular matrices ip' — UiipU^. Indeed, suppose it is possible and set ip = 1. Then U1U2 = e la I, which 
obviously is not true if a ^ 0. 
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The integral over A can be taken in an explicit way, because we can (under the 
condition p > 0, which follows from the definition p = tr((p^(p)) close the contour in 
the upper half-plane, and inside the contour there's a pole A = | of order N 2 . The 
integral equals 2m x (residue in this pole). Thus, 

e lXp 1 2tt /p^ 2 



We obtain 



■'BV 



N 2 r 1 

7T / , 1 ^2 N 2 F{ P )_ £ 



{I p - 1)1 J <V" V < 2 - 6 > 

Next we make a change p — > iV 2 p to reduce the integral to the case, in which the 
steepest descent method is applicable: 



{N2n)N2 1 dp-e N2 ^^n P )~i). (2.7) 



* BV (iV 2 -l)! J ^ p 
The stationary point equation is 

Po: --~ + F'(p o ) = 0. (2.8) 
po 2 

Now we specify the potential: F(p) = — fp 2 - Then the equation takes the form 

- - ~ - § Po = 0. (2.9) 
Po 2 2 

A positive solution of the equation is 

Po = ^( V / lT%-l) (2.10) 
We get the following expression for the partition function in the large N limit: 

Zbv= ^tl{ * yy<^-?-4> (2 .n) 



(iV-l)!^^ 9 (l + ri _) 



Using Stirling's formula for the factorial (n! ~ (- ) n \ / 2im), 

Z BV = (e ■ vr)^! f f + -1—) ^pd-*-,* ( 2 .12) 
po V 2 2-po/ 

From (|2.10| ) it is clear that in the limit g — > oo : p ~ (^) x ^ 2 ^ 0. Therefore the 
only growing contribution to the free energy is given by the logarithm, and 

S -> iln(^). (2.13) 

9->oo 2 
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3. Solution of the D=0 matrix model with tr((<^<y2) 2 )-interaction 



In this section we consider a model which has the N x N complex matrix tp as its 
dynamical variable. We will consider the partition function 



Z = J d<p eX p[--tr(^V) - ^ tr ((<^V) 2 )]- 



(3.1) 



Correlation functions of operators tr((<^ ip) m ) were obtained in [4]. All planar Green's 
functions are constructed with the help of the eigenvalue distribution function u(X), 
defined on a segment A G [0; y/z\ (here A should be interpreted as the modulus of an 
eigenvalue, so A > 0), which has the following form: 



u(X) 



A 2 



(3.2) 



where z = ^(\/l + 24g — 1). Note that as g — > oo z — > 0. At the same time the area 
under the plot u(X) is always equal to unity: J u(X) dX = 1. As a result, obviously 



lim u(X) = 5(A), 



(3.3) 



where the limit should be understood, of course, in the sense of distributions. One 
can see that in the limit g — > oo all eigenvalues tend to approach zero. Never- 
theless, it's not possible to use directly the asymptotics ( |3 .3| ) in order to obtain 
the partition function at large g, as in the expression for the free energy there's a 
term J dXu(X) In A, and the logarithm is singular at A = 0. We propose another 
method for the solution of a complex matrix model, which is free from the men- 
tioned drawback, and, using this method, we will prove at D = the statement 
presented above about the equality of free energies. First we make in the integral 
the Hubbard-Stratonovich transformation of the second term in the exponent: 



Z 



d^ dip drj exp 



--tr(</?V) - tir] 2 + iJ -^tr(r^V) 



(3.4) 



where rj is a hermitian matrix. Now we make the gaussian integration over <p: 



Z 



drj exp 



— tr rj — Ntr In 




(3.5) 



After passing to the integration over (real) eigenvalues A& of a hermitian matrix rj 
the integral can be calculated by the stationary phase method. Upon varying the 
exponent over A& we get the following equations: 



A, 



2X k + 



jjX k 



k = 1...N. 



(3.6) 



; , ., -1/2 + u 

In spite of the fact that we're integrating over real space ~R N , the values A& in the 
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stationary point in our case are complex — it is a usual situation for the steepest 
descent method. Let's now pass to the continuous case (N — ► oo), changing A& — > 
vN\(x) (x G [0, 1] is an analog of the index k for the continuous case). Then we get 
a system 

— = A + ^ 3.7 

A-v 2 ^gX + \ 1 ; 

u(z/)di/ = 1, (3.8) 



c 



c 



where v(X) = |f. From the discrete form ( |3.6| ) of our equation one can see that, 
if {Xk}k=i...N is a solution, so is {—\* k }k=i...N, therefore the contour C is symmetric 
with respect to the axis Oy = ImA. Next we multiply the first equation by 1 — 2i v /gA 
and introduce a new function v(X) = (1 — 2i-JgX)v(X). According to the position of 
the contour, f c v(X) dX = D is real (because trA = f c Xu(X)dX is imaginary). Then 
we get an equation 

p J = + A(1 _ 2 .^_ A) (3 g) 

c 

We will look for a "resolvent" $(A) = f X C in the following form 

c 



$(A) = + A(l - 2i^A) + (AA + B) a/(A + 6*) (A - 6) (3.10) 

In the limit A — > oo the asymptotics is known $(A) — > y (which is clear from the 
definition of the resolvent and the normalization condition — see above), but since D 
is unknown we will use the following conditions: the coefficients of A 2 , A 1 , A vanish 
and J c v(u)du = J c ^i^u = ^ ^he nrs ^ three conditions can be rewritten in the 
following form: 



A 2 
A 1 
A 



-2i^/g + A = (3.11) 
l + B + i</g(b* -b) = (3.12) 
v^(Re6) 2 + B1mb- y/g = (3.13) 



In writing out these equations we took the positive value of the square root on 
the continuation of the segment [—b*; b] to the right of b. Then on the axis, which 
is perpendicular to this segment and intersects it in the middle, the square root 
is imaginary with a positive (above the intersection point) or negative (below the 
intersection point) imaginary part. Suppose Ao G C. Taking into account what has 
been said and the Sokhotsky-Plemelj formulas, 

lim ($(A + *e)-$(A -2e)) = 2i(A\ + B) y/ (\ Q + 6*) (6 - A ) = -2wnJ(A ), (3.14) 
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thus 



v(X) = --(AX + B)y/(X + b*)(b - A), A e C. 

71 



Therefore we can rewrite the last condition in the form 



1 

TT 



-6* 



Av + B 
1 - 2iy/gi/ 



y/(v + b*)(b-v) dv = 1. 



So we're interested in the integrals of the following type: 



(3.15) 



(3.16) 



b 

h = J yj(v + 6*) (6 - v) dv 



-b* 



o 

-St 



dv 



2iyfgv 



y/{v + b*){b-v). 



(3.17) 



(3.18) 



They can be calculated in a similar fashion. Actually these integrals are contour 
integrals along the upper bank of the cut from left to right. Since on the lower bank 
of the cut the integrand takes on an opposite value, the value of our integral is half 
the value of the integral along a contour closed clockwise. The integral along a closed 
contour in the first case I\ is just 2m x (residue at infinity), and in the second case 
I 2 one also has to take into account the residue at the pole. Adopting the notation 
r(iz) = ^(j, + &•)(&-„), 

h 

I 2 = 7i (res 



7TTeS{T{V), OO J 

t(v) 



1 - 2iy/gv 



, oo) + res( 



r{v) 



1 - 2iy/gv 



,^o)), 



(3.19) 
(3.20) 



were v 



-^q= is the position of the pole. Upon calculating we obtain: 



|(Re6) 2 



h = 



7i 



2Va 



(Reb) 2 + (Im6 + 



?V9 



1/2 



Im6 



2^ 



(3.21) 



(3.22) 



Note that from the definition of the integral hig) it is clear that it is nonsingular 
as g — > 0. Formally there are singularities in the explicit expression obtained for J 2 , 
but it is easy to see that they cancel each other. So it is a check of our calculations. 
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Let's now write out all the conditions on (B, Re b, Imb), eliminating A: 



B = -l-2y/gTmb (3.23) 

= y/g~(Reb) 2 + B Imb - y[g~ (3.24) 

1 2,/glmb 

1 = -h+ J 2 (3.25) 

7T 7T 

The equations are real, hereupon the number of unknowns coincides with the number 
of equations. This confirms that it was quite relevant to take into account the 
symmetry of the contour. Next we introduce the parameter e = ^= and eliminate 
B: 

1 = ^(Re&) 2 - (Im&) 2 - elmb (3.26) 
1 = i(Re&) 2 - (Imb) 2 -elmb + Imb \/(Reb) 2 + (Imb + e) 2 . (3.27) 

Passing to the limit e — > 0, introducing b = \imb(e) and eliminating (Refe) 2 with the 
help of the first equation, we get 

~ = Imb ^3(lmb) 2 + 1. (3.28) 



The solution is as follows: 

Im6=-J=; Re6=-^=. (3.29) 

Thus, the contour C lies above the real axis, so the pole z/ = —is never hits the 
contour, as one could expect. The ends of the contour in the limit e — > turned out 
to be on a finite distance from each other as well as from the origin. It means that this 
approach correctly reflects the asymptotic properties of the model — "collapsing" 
of the contour, which took place in (|3.3j ) as g — ► oo, is absent. Let's now return to 
the problem of calculating the partition function. As it is clear from (|3.5|) , the free 
energy in the limit iV —>■ oo is: 

S = j d\v{\) ^A 2 + ln(^ - iy^gX^j - J d\dnv{\)v{n) In |A - n\. (3.30) 

From the presented results one can see that there's a finite limit of the function v(X), 
as g — > oo. Therefore, the main contribution to 8 in the described limit is given by 
y/g, which appears under the first logarithm, thus 

£ - -In*?, (3.31) 

g—foo Z 



which exactly coincides with the analogous expression for the vector model (|2.13| ) 
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4. Discussion and conclusion 

In this paper we have presented a detailed consideration of the D = complex 
model. In particular, we obtained an expression for the partition function by a 
method, analogous to the method presented in [3]. To the leading order in g in 
the strong coupling limit the hypothesis about the analogy of the matrix and vector 
models turns out to be true. Moreover, it is also true in the next-to-leading order if 
one uses in the vector model a finitely renormalized charge g = a ■ g (a = const.) 
instead of g. This is clear from the general structure of the free energy for both 
models: 

E= l -\ng + (3 Q + ^ + ... (4.1) 

We can always set the constant (3q to any value by changing in an appropriate way 
g — > a ■ g. What is nontrivial is the fact that the functional structures in the 
strong coupling expansion of the matrix and vector models coincide, although the 
coefficients may differ (a few first coefficients can be made equal by an appropriate 
finite renormalization of parameters). 

Our hypothesis is also true at D — 1, which can be seen from the explicit form 
of the free energy asymptotics for the matrix and vector models: both behave like 
E ~ k ■ g 1 / 3 (k = const). The question of whether our hypothesis is true for D > 2 
remains open. 
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